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I.  INTRODUCTION 

The  main  objective  of  this  work  was  to  determine  02(*A)  yield 
from  the  nuclear  pumping  of  a  noble  gas-02  mi‘xture.  To  achieve  this  goal 
it  was  necessary  to  benchmark  a  chemical  kinetics  model  to  predict  the 
scaling  of  the  02(k)  yield  with  noble  gas  pressure  and  with  thermal 
neutron  flux.  A  second  goal  was  to  test  the  suitability  of  nuclear 
produced  02(*A)  for  the  operation  of  an  iodine  laser.  A  third 
objective  was  to  explore  the  potential  of  other  methods  for  nuclear 
generation  of  02(*A).  This  led  to  the  Identification  of  photolytic 
dissociation  of  ozone  as  a  promising  alternative  to  direct  pumping  of 
02(1*). 

A  brief  description  of  nuclear  pumping  is  illustrated  by  Fig.  IA.  The 
system  typical  of  this  set  of  experiments  was  an  aluminum  tube  internally 
coated  with  Boron-10.  The  Boron-10  intercepts  thermal  neutrons  escaping 
from  the  TRIGA  reactor  core  and  undergoes  a  nuclear  reaction  to  generate 
MeV  alpha  particles  and  lithium  ions.  These  ions  stream  through  the  gas 
creating  ionization  and  excitation  as  they  slow  down,  resulting  in  the 
formation  of  a  plasma.  An  alternative  approach  to  pumping  the  gas  is  shown 
in  Fig.  l.B  where  the  medltui  contains  a  gaseous  pumping  source  i.e.,: 

3He,  UFg  or  BF3.  In  the  flashlamp  studies  proposed  later  in  Section 
1 1.4,  the  photolytic  fluorescer  requires  a  volume  source  to  achieve  the 
required  Intensities.  3He  is  suggested  for  preliminary  experiments  since 
It  is  the  easiest  to  handle.  In  a  volume  source,  the  thermal  neutrons  are 
captured  by  the  reactive  medium,  resulting  in  nuclear  reactions  that 
generate  MeV  ions  directly  in  the  gas.  The  advantage  of  the  volume  source 
over  a  surface  source  is  that  virtually  all  of  the  energy  released  is 
deposited  In  the  gas  while  with  the  surface  source  only  a  fraction  of 
energy  Is  deposited  because  ions  travelling  in  the  wrong  direction  are 
absorbed  by  the  coating  or  the  tube  wall. 


Laser  medium 


A.  Coated  tube  design 


B.  Volume  source  design 

Figure  1  Nuclear  pumping  sources  are  of  two  forms;  A)  a 
surface  source  where  the  high  energy  ions  are 
created  in  a  coating  on  the  walls  of  the  cell, 
or  B)  a  volume  source  where  a  gas  additive  is 
used  to  create  the  high  energy  ions  in  the 
plasma. 


In  the  following  section  we  first  discuss  the  results  that  lead  to 
the  proposal  for  02(*A)  production  by  nuclear  photolysis  of  ozone.  In 
Section  III  experimental  data  obtained  from  measurement  of  direct  nuclear 
pumping  of  (^(^A)  are  presented.  This  includes  results  from  transport 
measurements,  parametric  studies,  and  kinetic  modeling.  Also  a  summary  of 
operating  parameters  observed  for  a  nuclear  sustained  electrical  discharge 
is  included.  This  device  was  developed  for  steady-state  generation  of 
02(*A)  and  was  used  in  the  transport  studies. 

II.  o2(1a)  production  by  ozone  photolysis 

A  key  discovery  during  this  work  was  the  realization  that  nuclear 
pumping  favors  the  generation  of  ozone,  particularly  at  higher  pressures. 
In  order  to  capitalize  on  this  result,  we  considered  the  photolytic 
decomposition  of  ozone  into  02(^a)  as  an  alternative.  The  peak  ozone 
yield  observed  in  a  preliminary  experiment  was  approximately  0.5  Torr  O3 
in  a  4000  Torr  He,  20  Torr  O2  mixture  at  the  observation  cell.  Based  on 
a  computer  model  for  the  decomposition  of  ozone  by  a  KrF*  fluorescor  (x  = 
249  nm),  it  appears  to  be  possible  to  obtain  >  5  Torr  02(*A)  starting 
with  the  0.5  Torr  of  O3  cited  above.  (Taking  into  consideration  the  gas 
expansion,  then  between  1  to  2  Torr  peak  O3  exists  in  the  irradiation 
zone.)  The  computer  model  indicates  that  with  an  initial  O3  pressure  of 
2  Torr,  ~7  Torr  of  02(*A)  can  be  generated,  representing  a  35% 
concentration  with  respect  to  the  total  oxygen  pressure.  Furthermore,  it 
appears  that  the  nucl ear-pumped  flashlamp  intensities  necessary  to 
completely  dissociate  the  ozone  are  within  current  technology. 

Figure  2  shows  a  straight  flow-through  design  which  would  use  a  rela¬ 
tively  fast  flow  rate  to  generate  the  ozone  in  the  first  half  of  the 
reactor  pulse  and  photolyze  the  O3  into  02^A)  during  the  second  half 
of  the  pulse.  Based  on  preliminary  experiments,  the  nuclear  production  of 
ozone  appears  to  be  relatively  efficient.  The  main  production  mechanism 
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that  appears  to  be  responsible  for  ozone  production  is  three-body 
recombination  of  ground  state  atomic  oxygen  with  ground  state  molecular 
oxygen  In  the  presence  of  a  third  body  such  as  helium  or  oxygen. 

1.  Ozone  Measurements 

The  experimental  apparatus  that  was  used  to  determine  the  ozone 
concentration  is  shown  in  Figure  3.  The  reaction  zone  consisted  of  a 
1.27  cm  diameter  Boron-10  coated  aluminum  tube  60  cm  long  which  was  placed 
adjacent  to  the,  TRISA  core.  The  excited  He-Oz  mixture  was  then  flowed 
through  10  meters  of  tubing  (polyflow)  to  an  absorption  cell  where  the 
ozone  concentration  was  measured.  The  detection  of  O3  was  accomplished 

O 

with  a  He-Hg  lamp  radiating  at  2537  A,  monitored  by  an  RCA-31034 
photomultiplier  coupled  to  a  Heath  monochromator.  In  the  case  of  the  high 
pressure  studies,  the  gas  was  simply  vented  out  of  the  tube.  At  pressures 
below  1  atmosphere,  an  oxygen  serviced  vacuum  pump  was  employed. 

Experimental  results  from  this  system  are  shown  in  Figures  4  and  5. 

In  Figure  4  the  ozone  pressure  increases  relatively  rapidly  with  increasing 
He  pressure  and  appears  to  saturate  at  around  3500  Torr  of  He.  This 

saturation  is  attributed  to  power  deposition  saturation.  This  problem  can 
be  overcome  with  a  smaller  diameter  tube  that  allows  the  alpha  particles  to 
reach  the  centerline  at  desired  maximum  pressure.  However,  this  was  not 
attempted  in  present  experiments.  The  solid  line  in  this  figure  shows  the 
actual  experimental  data  (approximately  0.5  Torr  of  ozone  at  the 
observation  point  10  meters  from  the  initial  excitation  zone).  If  the 
expansion  of  the  excited  gas  slug  is  taken  into  account,  the  equivalent 
ozone  pressure  in  the  excitation  zone  in  the  core  is  estimated  to  be 


Experimental  apparatus  used  to  monitor  0-  yields  from  nuclear  pumping. 


HELIUM  PRESSURE  (Tori 

Experimental  results  obtained  with  the  apparatus  diagram* 
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approximately  2  Torr.  These  0^  yields  are  more  than  more  than 
sufficient  for  the  (^(^A)  generation  technique  suggested  in  the  follow¬ 
ing  section.  Results  shown  in  Figure  5  were  obtained  for  cw  operation  of 
the  reactor,  eliminating  the  uncertainty  associated  with  expansion  of  the 
gas  slug.  The  increase  of  ozone  yield  with  power  deposition  is  clearly 
demonstrated  by  these  results. 

2.  Ozone  Photolysis  Calculations 

In  the  photolytic  generation  of  O2 ( ^ from  ozone,  a  rather  unique 
process  can  be  exploited.  First,  the  O3  is  rapidly  decomposed  by  UV 
light  (typically  less  than  300  nanometers  and  greater  than  230 
nanometers)  into  02(*A)  and  0(^0).  The  0(^D)  is  rapidly  quenched 
by  other  species  to  ground  state  atomic  oxygen  which  can  then  undergo 
three-body  recombination  with  ground  state  molecular  oxygen  and  a  third 
body  (e.g.  He)  to  form  ozone  (see  Figure  6).  This  mechanism  forms  a 
positive  feedback  which  enhances  the  02^)  concentrations  by  as  much 
as  an  order  of  magnitude  compared  to  the  initial  ozone  concentrations. 

This  technique  operates  best  at  high  pressures  where  the  three-body 
recombination  process  is  fast  throughout  the  pulse. 

Note  that,  as  shown  in  Figure  6,  the  ozone  is  initially  photolyzed 
into  (^(^A)  and  0(*D)  or  0(^p)  and  02(%).  In  the  present  calculation  a  90% 
branching  ratio  was  assumed  for  the  02(*A)  -  0(^0)channel .  This  is 
a  rather  conservative  estimate  since  much  of  the  recent  literature  states 
that  this  branching  ratio  is  approximately  unity  (see  Appendix  A). 

A  fairly  detailed  kinetic  model  has  been  developed  for  calculation  of 
O3  photolysis  yield.  This  model  is  described  in  more  detail  in  Appendix 

Results  from  these  calculations,  shown  in  Figure  10,  indicate  that  as 
the  flashlamp  power  increases,  the  ozone  concentration  slowly  decreases. 


i.  a! 1m  * 


1IZ  Torr  O3,  18  Torr  Og  ' 

3800  Torr  He 

PEAK  FLASHLAMP  INTENSITY 
=  600  W/cm2 
PULSE  WIDTH  =  13.8  ms 

__j _ 1 _ 1 _ 1 _ 

0  0.04  0.08  0.12  0.16  0.20 

TIME  (seconds) 

Fig.  7  Time  dependent  solution  of  chemical  kinetics 
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The  peak  flashlamp  intensity  assumed  was  600  W/cm2  which  appears  to  be 
achievable  with  current  nuclear  pumping  technology.  The  rise  in  the  atomic 
oxygen  concentration  is  a  direct  result  of  the  photolytic  decomposition  of 
the  ozone.  The  (^(^A)  concentration  exceeds  the  ozone  concentration 
almost  immediately  due  to  the  positive  feedback  mechanism  already  outlined. 
This  calculation  suggests  that  starting  with  2  Torr  of  ozone,  up  to  7  Torr 
of  (^(^A)  can  be  produced.  This  represents  an  (^(^a)  concentration 
equal  to  35%  of  the  total  oxygen  concentration,  well  above  the  requirements 
for  pumping  an  iodine  laser. 

The  next  several  figures  show  the  O2 ( ^ A)  yield  as  a  function  of 
partial  and  total  mixture  pressure,  flashlamp  intensity  and  duration.  This 
parametric  study  indicates  that  a  flashlamp  peak  intensity  of  150  W/cm2 
is  all  that  is  necessary  to  create  the  (^(^a)  concentration  (>  17%) 
required  for  operation  of  an  iodine  laser. 

The  most  interesting  point  of  Figure  8  is  that  if  the  flashlamp 
intensity  falls  between  600-1000  W/cm2,  this  technique  can  generate 
approximately  7-8  Torr  of  02( ^ A)  from  a  mixture  of  18  Torr  O2.  2  Torr 
O3,  and  3800  Torr  He.  For  a  reduced  ozone  pressure  (0.5  Torr  ozone)  and 
the  above  flashlamp  intensity,  between  5  and  6  Torr  of  (^(^A)  will  be 
generated.  Thus,  there  is  a  broad  range  over  which  significant  (^(^A) 
yields  can  be  obtained. 

3.  System  Parametric  Studies 

We  have  found  that  the  three-body  recombination  rate  for  ozone 
formation  increases  directly  with  He  pressures.  Thus  the  enhancement 
mechanism  for  02(*a)  production  from  the  photolytic  decomposition  of 
ozone  increases.  This  is  apparent  in  Figure  9  which  graphs  the  (^(^A) 
yield  as  a  function  of  the  He  pressure  and  as  a  function  of  various 
flashlamp  Intensities.  It  is  important  to  remember  that  as  the  flashlamp 
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Fin.  8  Parametric  study  of  the  02 ( ^ A )  yield  as  a  function  of 
the  flash! amp  intensity  for  various  O^tOg  mixtures. 
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intensity  increases  the  depletion  of  the  ozone  during  the  pulse  increases. 
Since  ozone  strongly  quenches  02^a)»  the  0£( ^A)  lifetime  Increases. 

As  a  result,  higher  yields  of  02(*a)  are  achieved  with  higher  flashlamp 
intensities.  This  point  is  also  brought  out  in  Figure  10  where  with  a  1 
KVI/cm^  peak  flashlamp  intensity,  7.5  Torr  (^(^A)  with  an  [(^(^A)]/^^!)] 
ratio  of  .58  is  generated.  The  effective  half life  of  the  generated 
(^(^A)  is  70  ms.  A  further  increase  in  the  effective  lifetime  of 
(^(^A)  is  possible  using  a  more  intense  UV  source. 

Figure  11  shows  the  dependence  of  02( ^A)  yield  on  flashlamp 
duration  for  different  02^3  mixtures  and  for  a  constant  peak  flashlamp 
intensity.  The  production  efficiency  of  (^(^A)  is  large  for  shorter 
pulses  because  of  the  rapid  decomposition  of  O3  which,  as  already  noted, 
is  a  strong  quencher  of  02(^A). 

Thus,  from  this  parametric  study  we  conclude  that  02( ^ A)  production 
by  this  method  is  suitable  for  iodine  laser  operation.  Experimental 
testing  of  this  model,  however,  is  still  needed  in  order  to  verify  these 
results  and  determine  optimum  conditions. 

4.  Nucl ear-Pumped  Flashlamp  Concepts 

A  nucl ear-pumped  flashlamp  concept  for  use  in  this  scheme  is 
Illustrated  in  Figure  12.  In  this  approach  the  high  power  densities  needed 
at  the  center  tube  section  (where  the  flowing  ozone-helium  mixture  passes) 
are  achieved  by  geometric  focussing  using  a  highly  polished  aluminum 

surface.  This  focussing  concept  concentrates  the  light  generated  in  the 
large  fluorescor  volume  while  transferring  it  to  the  smaller  volume  of  a 
very  highly  absorbing  medium,  in  this  case  ozone.  A  possible  fluorescor  is 
KrF*  which  emits  at  a  wavelength  of  249  nanometers,  roughly  the  peak  of  the 


PARTIAL  PRESSURE  (Torr) 


FOCAL  POINT 


Fig.  12  A  nuclear  pumped  flashlamp  design  which  can  focus  the 
fluorescence  from  a  large  volume  to  a  smaller  volume 
resulting  In  a  net  Increase  In  the  intensity  in  the 
center  cell. 
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ozone  absorption  cross  section  (~  11  megabarns).  The  nuclear-induced 
fluoresence  efficiency  defined  as  ratio  fluorescence  energy  to  energy 
deposited  in  the  fluorescence  media  has  been  calculated  to  be  approximately 
14%  for  KrF*,  a  reasonably  high  value.  The  flashlamp  intensity  is  a  function 
of  the  power  deposited  in  the  fluorescor,  the  fluorescor  efficiency,  the 
coupling  efficiency  of  the  fluoresence  to  the  inner-cell  and  the  ratio  of  the 
volume  of  the  fluorescor  to  the  surface  area  of  the  ozone  cell.  Thus: 


rn  ... _  r„  -i  [Volume  of  Fluorescor] 

[Pin  (fluorescor)]  COupl i ng^  [Surface  area  of  absort 


Then,  for  a  TRIGA  reactor  which  has  a  thermal  neutron  flux  of  approximately 
2.5x10*5  n/cm2s,  the  power  deposition  in  5  atms  of  3He  will  be 
approximately  120  W/cm3.  If  a  parabolic  reflector  is  assumed  with  a 


coupling  efficiency  of  80%  to  the  inner  cell,  then  for  a  flashlamp  intensity 
of  1  KW/cm2,  the  ratio  required  for  the  flashlamp  volume  to  the  surface 


area  of  the  ozone  cell  is  -75. 


5.  Conclusions 

In  conclusion,  the  nuclear  generation  of  02(*a)  by  the  photolysis 
of  ozone  appears  very  promising.  Indeed,  a  TRIGA  reactor  neutron  pulse  can 
generate  workable  yields.  The  theoretical  scheme  examined  used  a  KrF-excimer 
flashlamp  to  photodissociate  the  ozone.  This  generates  OgC ^A)  directly 
along  with  atomic  oxygen.  The  latter  results  in  a  feedback  mechanism  which 
greatly  enhances  Og^A)  yields.  For  the  case  of  0.5  Torr  O3,  19.5 
Torr  02»  and  3800  Torr  He  (corresponding  to  experimental  ozone  yields)  and 
a  1  kW/cm2  flashlamp,  the  02(*a)  yield  will  be  approximately  10  Torr, 
corresponding  to  30%  of  the  total  O2  concentration.  The  associated 
half-life  of  the  02(*A)  is  estimated  to  be~70ms. 


The  idea  which  has  just  been  discussed  is  a  direct  result  of  the 
experiments  that  were  conducted  on  the  direct  nuclear  pumping  of  Ar-02  to 
form  O2 ( ) -  These  studies  revealed  that  the  most  important  species 
present  in  the  afterglow  is  ozone.  Consequently,  the  scheme  just  described 
would  take  advantage  of  the  large  ozone  yield.  The  work  which  will  now  be 
discussed  involved  initial  experiments  on  direct  pumping  of  (^(^A)  and 
the  associated  development  of  reliable  diagnostics  for  these  and  subsequent 


experiments. 

The  techniques  developed  used  a  flowing  afterglow  design  with  a 
Boron-10  coated  excitation  cell  coupled  to  a  6-meter  flow  path  of  pyrex 
tubing.  The  tubing  transported  the  excited  gas  slug  from  the  radiation 
zone  to  the  detectors  which  were  located  behind  concrete  shielding  to 
reduce  radiation  noise.  Another  diagnostic  developed  was  an  in-core  fiber 
optic  designed  to  be  highly  resistant  to  both  radiation  damage  and 
radiation-induced  fluorescence.  These  techniques  were  then  used  to 
determine  both  03  (*A)  and  O3  yields  from  a  direct  nuclear  pumped 
system  and  also  from  a  hybrid  nuclear -electrical  system.  Emissions  at  1.27 
microns  and  at  6340  A  were  used  to  determine  the  O2 ( )  concentration. 

O 

The  ozone  was  simultaneously  monitored  by  the  absorption  of  the  2537-A  line 
from  a  mercury  discharge  directed  through  an  observation  cell  approximately 
108-cm  long.  These  results  were  applied  as  the  benchmark  to  the  modeling 
which  was  then  used  to  estimate  how  the  02(k)  yield  scales  with 
increasing  noble  gas  pressure  and  also  with  the  reactor  power  levels. 


1.  Experimental  Arrangement 


The  experimental  apparatus  is  diagrammed  in  Figure  13.  The  O2 ( ^ 


was  generated  in  an  aluminum  tube  coated  on  the  inside  with  Boron-10  and 


placed  adjacent  to  the  TRIGA  reactor  core.  The  argon  and  oxygen  mixt'-  js 


IRRADIATION 
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Fig.  13  Experimental  set-up  for  monitoring  CL  (  A)  at  1.27  m.  Also, 
the  0-  density  is  determined  with  the  absorption  of  2537  A 
light Jfrom  a  He-Hg  lamp. 


were  flowed  through  this  excitation  cell  into  a  6-meter  long  pyrex  vacuum¬ 
line  which  transported  the  gas  from  the  radiation  zone  to  the  detectors. 

The  flow  velocities,  on  the  order  of  20  meters  per  second,  were  provided  by 
a  roots-type  pump.  A  LN2  cooled  germanium  detector  with  an  interference 
filter  centered  at  1.27  microns  was  used  to  observe  the  emissions  from 
02(*A).  A  second  detector,  also  a  germanium  unit,  was  used  without  a 
filter  to  monitor  broadband  emissions.  By  monitoring  the  emissions  from 
900  nm  -  1.4  pm  and  comparing  their  behavior  to  the  1.27  pm  line  it  is 
possible  to  determine  how  much  the  NO  +  0  afterglow  (caused  by  N2 
impurities)  affects  the  O2 ( )  measurements.  Also  depicted  in  this 
figure  is  the  helium  mercury  lamp  which  illuminated  the  observation  cell 
with  UV  light  (2537&).  Ozone  concentrations  were  determined  by  absorption 
of  this  light  as  monitored  by  a  Heath  monochromator  coupled  to  an  RCA  31034 
photomultiplier.  The  shielding  shown  was  necessary  to  prevent  radiation 
from  the  reactor  pulse  from  reaching  the  germanium  detectors,  which  are 
highly  susceptible  to  interference  from  gamma  and  x-ray  radiation. 

2.  Detector  Calibration 

The  detector  calibration  technique  used  a  relatively  straightforward 
method  that  involves  a  flowing  afterglow  with  an  NO2  titration.  The 
technique  relied  upon  the  following  reactions: 

1)  NO2  +  0  -*■  NO  +  O2  (very  fast) 

2)  NO  +  0  N02  +  hv  (slow) 

and  the  derived  relation: 


3)  kdet  =  K  MM. 


NO 


where  K^y  =  detector  calibration  constant 

K  =  literature  constants 
[0]  =  atomic  oxygen  concentration 
[NO]  =  nitric  oxide  concentration 

'no  ■ 


observed  intensity  of  NO-O  glow. 


The  above  set  of  reactions  is  rather  unique  because  the  large  difference  in 
their  reaction  rates  allows  for  a  simple  method  of  determining  the  atomic 
oxygen  concentration.  Suppose,  for  example,  the  NO2  and  atomic  oxygen 
concentrations  are  made  equal.  Reaction  1  will  very  quickly  exhaust  the 
supply  of  atomic  oxygen.  Since  there  is  no  regenerative  mechanism  for 
atomic  oxygen  in  this  system  once  it  is  gone,  both  reactions  will  cease. 

As  a  result,  the  light  emitted  by  reaction  2  will  "disappear."  So  by 
reversing  the  above  argument,  the  termination  of  the  NO  +  0  afterglow 
represents  the  titration  end  point  where  the  NO2  and  atomic  oxygen 
concentrations  are  equal. 

The  experimental  apparatus  used  for  the  calibration  is  shown  in  Figure 
14.  As  the  O2  passes  through  the  microwave  cavity,  it  forms  atomic 
oxygen,  02(*A),  and  various  other  species.  The  flow  tube  for  viewing 
the  afterglow  is  isolated  from  the  discharge  by  a  light  trap  downstream  of 
the  microwave  cavity.  The  NO2  is  injected  with  a  small  port  upstream  of 
the  infrared  detector  which  views  the  1.27-micron  emissions  from  both  the 
NO2  glow  and  the  (^(^A).  The  contribution  from  0 2 C ^ A)  is  subtracted  out 
of  the  signal  through  an  initial  background  measurement  (no  NO2).  This 
allows  a  direct  measurement  of  INg.  Then,  by  performing  an  NO2 
titration  as  described  above,  the  atomic  oxygen  concentration  is  readily 
obtained.  Finally,  since  reaction  1  is  so  fast  compared  to  2,  the  NO 
concentration  can  be  taken  as  the  initial  NO2  concentration. 

The  NO2  flow  rate  is  then  measured  using  a  calibrated  volume 
technique.  At  the  appropriate  NO2  flow  rate,  the  flow  is  diverted  into 
an  evacuated  calibrated  volume.  The  time  required  for  the  calibrated 
volume  to  fill  to  a  given  pressure  gives  the  flow  rate  of  the  NO2.  This 
technique  is  very  reliable  since  the  relatively  long  fill  times  involved 
result  in  a  high  accuracy  in  the  flow  rate  determination. 


Figure  15  illustrates  concentration  changes  during  operation  of  this 
system.  A  computer  simulation  of  the  atomic  oxygen  decay  was  performed 
and  compared  to  experimental  results  to  obtain  this  graph.  Once  the  NO2 
was  injected,  it  was  converted  directly  into  NO  and  the  atomic  oxygen 
density  was  quenched  abruptly  as  it  was  consumed  by  the  NO2  +  0  reaction. 
Using  the  computer  model  to  estimate  this  jump,  it  was  possible  to  predict 
the  atomic  oxygen  and  NO  concentrations  15  cm  downstream  of  the  NO2 
injector.  These  results  were  then  used  in  Eq.  (3)  to  obtain  K^y  and 
listed  in  Table  1.  These  data  have  been  reproduced  under  several  different 
oxygen  pressures. 

3.  Transport  Measurements 

To  minimize  the  effects  the  radiation  has  on  the  detectors  used  it  was 
apparent  that  the  transport  studies  can  only  be  completed  at  low  steady 
state  reactor  power  levels.  However,  at  reactor  powers  of  1  MW  or  less, 
the  02(1A)  concentration  was  considerably  less  than  what  was  obtained 
from  a  pulse.  Consequently,  a  nuclear  radiation  sustained  discharge 
(Figure  16)  was  constructed  to  enable  steady  state  production  of  the  O2 ( ) 
densities  observed  during  a  pulse.  By  utilizing  the  nuclear  reactor  as  an 
external  ionization  source  0 2 ( ^A )  densities  of  up  to  1.4  mTorr  in  1  Torr 
of  O2  were  generated.  Since,  the  radiation  levels  were  much  lower 
during  these  experiments,  the  detectors  could  be  moved  along  the  flow  path 
to  obtain  the  decay  rate  of  the  0 2 ( ^A ) .  By  matching  O2 ( )  and 
O3  yields  from  this  system  with  the  concentrations  typical  of  the  direct 
nuclear  pumped  approach,  it  was  possible  to  determine  the  (^(^A)  decay 
rate  for  the  conditions  observed  in  an  Ar-02  nuclear  pumped  system. 

The  results  of  the  transport  study  are  shown  in  Figure  17.  For  the 
mixture  of  150  Torr  Ar,  2  Torr  O2  and  1  millitorr  O3,  the  O2 ( ) 
lifetime  was  found  to  be  approximately  0.3  seconds.  The  lifetime  does  not 


injector  detector 


15.  Computer  simulation  of  reactions  occurring  after  mixing 
NOg  with  an  atomic  oxygen  stream. 
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Table  1 


Detector  Calibration  Constants  Obtained 
with  NO2  Titration  Experiments 


ADC  MODEL  403  @  1.27  m  [O^aJ  -  (1.5+.5)  x  1010 
RCA  31034  9  1.8kV  @  6340A  [021a3  -  (7.8+.5)  x  1010 


/molecules! 


\ 


cm3-mV 


/molecules\ 


Table  2 


02(1A)  Yields  from  a  Direct  Nuclear- Pumped 
'  Noble  Gas  -02  Mixture 


150  Torr-Ne,  2  Torr-02 
150  Torr-Ar,  2  Torr-02 
150  Torr-He,  2  Torr-02 


0.18  mTorr  02(^A) 
0.15  mTorr  02(^A) 
0.14  mTorr  02(^A) 
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10 B  COATED  PLATES 


Fig.  16  Hybrid  Nuclear/Electrical  Discharge  apparatus 
A  10B  coating  was  placed  on  an  aluminum 
plate  and  suspended  by  the  Plexiglas  walls. 


vary  significantly  for  the  three  cases  shown  because  the  dimole  reaction, 
(^(^A)  +  02^A)  -*■  02(^A)  +  O2,  was  not  important  at  the  low 
02(^A)  concentrations  involved. 

Figure  18  shows  the  measured  spread  of  the  excited  gas  slug  as  it 
passed  the  detectors.  A  gaussian  shaped  slug,  approximately  two-feet  long 
was  expected.  8ut  as  seen  in  these  diagrams,  for  the  case  of  helium  the 
excited  gas  slug  was  considerably  broader  than  the  initial  dimensions  of 
the  Boron-10  coated  aluminum  tubes  inside  the  reaction  zone.  Similar  re¬ 
sults  were  observed  in  Argon  and  Neon  but  the  magnitude  of  the  spreading 
was  less  pronounced.  The  differences  appear  to  be  largely  associated  with 
variations  in  gas  viscosity.  After  correction  for  this  expansion, 

02(*A)  measurements  indicated  the  neon  oxygen  mixtures  produced  the 
highest  yield.  These  yields  are  shown  in  Figure  19  (also  see  Table  2) 
where  the  neon-oxygen  mixtures  produced  0.18  millitorr  of  02(*A),  the 
argon-oxygen  mixtures  yielded  0.15  millitorr  02  (*A),  and  the  helium- 
oxygen  system  peaked  at  0.14  millitorr  of  O2  (*A).  These  results  were 
obtained  with  a  pressure  of  150  Torr  noble  gas  and  2  Torr  oxygen. 

4.  Parametric  Studies 

The  O2  (*A)  yields  were  observed  to  increase  with  the  buffer  gas 
pressure  and  with  increased  neutron  flux.  However,  as  stressed  earlier,  at 
high  buffer  gas  pressures,  ozone  production  is  favored.  Figure  20  shows 
the  scaling  of  the  02(*A)  yields  with  buffer  gas  pressure.  Notice  that 
in  argon  the  yields  appear  to  saturate  at  approximately  100  Torr  which  sug¬ 
gests  the  alpha  particles  have  reached  their  maximum  range  for  this  pres¬ 
sure.  In  the  case  of  neon-oxygen,  the  saturation  begins  at  a  considerably 
higher  pressure  (around  250  Torr).  This  is  to  be  expected  since  the  range 
of  alpha  particles  in  neon  is  larger  than  in  argon. 

The  UV  absorption  diagnostic  allowed  the  ozone  yields  to  be  monitored 
directly.  The  yields  were  observed  to  increase  with  O2  pressure,  thermal 
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Fig.  18  A  comparison  between  the  spread  out  excited  gas 
slug  and  the  absorption  path  used  to  monitor  0.,. 
The  expansion  factors  were,  He:3.2x,  Ne:1.5x 
Ar:l .2x. 
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Fig.  20  Scaling  of  O-^A)  yields  with  increasing  buffer 
gas  pressure;  He-02  mixtures  are  not  graphed 
due  to  the  low  02(la)  yields  obtained. 


neutron  flux,  and  noble  gas  pressure.  The  net  result  of  this  study  is  that 
helium  was  found  to  be  the  best  buffer  gas  for  producing  ozone  under  the 
present  conditions.  The  dissociative  resonant  energy  transfer:  He^  + 

O2  -*■  0  +  D+  +  2He  is  apparently  responsible  for  the  efficient  formation 
of  atomic  oxygen  in  the  discharge.  At  these  high  pressures  the  atomic 
oxygen  is  rapidly  converted  into  ozone  by  three  body  recombination. 

Figure  21  shows  the  uncorrected  yields  for  ozone  as  a  function  of 
increasing  O2  pressure.  The  apparent  decrease  in  ozone  yields  for  dif¬ 
ferent  noble  gases  was  merely  a  result  of  the  differences  in  spreading  of 
the  excited  gas  slug  as  it  was  transported  6  meters  to  the  detectors.  Upon 
correcting  for  this  expansion  (Figure  22),  the  helium  oxygen  mixture  is 
found  to  have  the  highest  yields.  The  scaling  of  the  ozone  yield  with 
noble  gas  pressure  is  diagrammed  in  Figure  23.  The  yields  in  the  He-02 
mixture  increased  at  a  faster  rate  than  in  Ne-Og*  which  in  turn  increased 
at  a  faster  rate  than  in  the  Ar-C^*  This  is  probably  due  to  the  more 
efficient  dissociative  resonant  energy  transfer  in  helium  which  results  in 
the  efficient  formation  of  atomic  oxygen. 

5.  In-Core  Fiber  Optic  Diagnostic 

The  in-core  diagnostic  system  which  was  developed  in  this  project 
employed  a  silica  lens  system  to  couple  an  observation  cell  with  a  fiber 
optic  bundle  (Figure  24).  The  silica  lens  system  was  necessary  to  locate 
the  fiber  optic  bundle  outside  of  the  extremely  intense  radiation  field 
which  occurs  close  to  the  reactor  core.  The  silica  lenses  were  chosen  to 
take  advantage  of  their  high  resistance  to  radiation  damage.  The  first 
series  of  measurements  were  made  with  a  high  purity  silica  fiberoptic 
bundle  which  had  an  attenuation  of  20,000  dB  per  kilometer  at  1.27  microns. 
Consequently,  it  was  not  possible  to  observe  emissions  from  02(*A)  at 


1.27  microns.  Since  the  attenuation  was  much  less  in  the  visible  regime 
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Fig.  23  Scaling  of  0,  yield  with  increasing  Noble  gas 
pressure  for  a  $3.00  TRIGA  pulse.  These  results 
are  also  corrected  for  the  expansion  depicted  in 
figure  17. 


(20dB/km),  measurements  were  possible  at  6340  A  but  this  required  elec¬ 
tronic  filters  to  suppress  the  fluorescence  from  the  Cerenkov  emissions 
generated  in  the  lens-fiberoptic  system.  However,  the  emissions  observed 
in  the  observation  cell  were  of  insufficient  intensity  to  produce  consis- 

o 

tently  reliable  results  at  6340  A.  To  overcome  this  problem  an  anhydrous 
silica  bundle,  which  has  an  attenuation  at  1.27  microns  of  only  20  db  per 
kilometer,  was  ordered.  This  would  have  allowed  operation  on  the  higher 
intensity  1.27  microns  emission,  however  there  was  insufficient  time  to 
perform  the  measurements  before  the  termination  of  this  contract. 

The  incore  apparatus  also  employed  a  flowing  afterglow  design  (Figure 
24).  With  this  design,  observation  of  the  afterglow  occurs  only  30  cm  from 
the  radiation  zone  vs.  6  meters  in  earlier  designs.  The  observation  cell 
was  Isolated  from  the  discharge  chamber  by  a  pair  of  light  traps  and  the 
light  transport  tube  (using  a  series  of  silica  lenses)  was  embedded  in 
borated  polyethelene  to  reduce  the  radiation  interacting  with  the  system. 

A  small  chopping  fan  was  used  at  the  entrance  to  the  light  transport  tube, 
i.e.  between  the  tube  and  the  observation  cell.  The  signal  exiting  the 
lens  fiber  optic  system  consisted  of  a  component  with  the  chopping  fre¬ 
quency  and  a  slowly  varying  emission  which  followed  the  neutron  pulse 
(w  -  100  Hz).  Therefore,  the  low  frequency  component  of  the  signal  was 
eliminated  with  a  100  Hz  high-pass  electronic  filter.  A  150  Hz  band  pass 
filter  was  then  used  to  eliminate  the  high  frequency  noise.  The  signal 
with  the  chopper  frequency  was  finally  reconverted  to  a  DC  level  with  a 
lock-in  amplifier  (PAR  H-8).  The  processed  signals  were  then  displayed  on 
the  Tektronix  scope  and  the  results  were  recorded  with  a  Polaroid  camera. 

6.  Kinetic  Models 

Modeling  of  the  argon  and  oxygen  nuclear  pump  system  has  advanced 
quite  significantly  since  the  beginning  of  this  project.  A  more  detailed 
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The  Light  Transport  Tube  allowed  relocation  of  the  fiber  bundle  outside  of  the  intense 
radiation  zone. 
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description  of  the  kinetic  model  and  some  of  the  calculations  performed 
with  it  are  presented  in  Appendix  B.  Calculations  with  the  model  exhibited 
the  same  dependences  with  increasing  noble  gas  and  O2  pressures  as  have 
been  observed  experimentally  (see  Figure  25).  One  of  the  more  significant 
points  in  the  modeling  is  that  the  results  indicate  that  the  O2 ( 
yields  will  increase  with  noble  gas  pressure  and  thermal  neutron  flux. 

As  the  O3  yields  begin  to  increase  with  the  noble  gas  pressure  and 
power  deposition,  the  O2 ( ^A)  lifetime  will  decrease.  This  is  expected 
since  ozone  and  02(*a)  mutually  quench  each  other  in  the  afterglow. 

Figure  25  shows  a  comparison  between  the  modeling  and  the  experiments  with 
increasing  noble  gas  pressure.  Relatively  good  agreement  is  generally 
observed.  However,  at  the  higher  argon  pressures,  the  ozone  concentration 
is  somewhat  overestimated  by  the  model.  This  may  be  because  of  a  power 
saturation  term  and  the  effect  of  oxidation  of  the  walls  have  not  been 
included.  (These  effects  are  now  being  added  to  the  model.)  The  results 
of  the  modeling  at  the  peak  of  the  pulse  are  important  since  they  indicate 
the  02(*A)  concentration  will  increase  with  argon  pressure  but  its 
lifetime  will  become  shorter.  Therefore,  less  O2 ( ^ A)  will  be  observed 
downstream.  ' 

7.  Nuclear  Sustained  Electrical  Discharge 

The  nuclear  sustained  discharge  was  constructed  to  perform  the  trans¬ 
port  measurements.  Throughout  these  experiments  this  system's  operation 
had  to  be  characterized  to  insure  that  the  conditions  present  in  a  direct 
nuclear-pumped  system  could  be  simulated.  The  apparatus  used  was  shown 
earlier  in  Figure  16.  It  was  constructed  out  of  a  plexiglas  box  in  which 
two  Boron-10  coated  aluminum  plates  were  supported  transverse  to  a  pair  of 
l/2-1nch  diameter  electrodes  by  the  outer  walls  of  the  Plexiglas  cell.  In 
operation  the  discharge  was  observed  to  be  unaffected  by  the  3oron  10 
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Fig.  25  The  lower  graph  shows  a  comparison  between  the  modeling 

and  the  experimental  results.  The  upper  graph  illustrates 
.the  equivalent  behavior  at  the  peak  of  the  pulse. 
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plates  which  are  isolated  from  both  the  anode  and  the  cathode  by  the 
Plexiglas  wails.  At  moderate  reactor  power  levels  (<1  MW)  a  relatively 
stable  discharge  was  obtained  which  could  operate  at  low  electric  fields 
(Figure  26).  We  have  observed  up  to  1.4  millitorr  of  0£(^a)  in  1  Torr 
O2  and  100-Torr  Ar  with  a  modest  electrical  power  deposition.  Figure  27 
demonstrates  how  scaling  to  higher  yields  is  indeed  possible  with  lower  gas 
pressures  and  higher  electrical  power  input.  This  graph  shows  that  the 
(^(^A)  concentrations  increase  relatively  rapidly  with  respect  to  the 
square  root  of  electrical  power  deposited  in  the  system.  If  the  curves  for 
150  Torr  argon  are  compared  to  100  Torr  argon,  the  slope  of  the  02^) 
yield  does  not  change.  However,  the  power  level  required  to  generate  a 
given  02^)  yield  is  reduced.  In  addition,  the  O3  yields  were 
observed  to  decrease  with  an  increase  of  the  electrical  power  input  (Figure 
28).  This  effect  is  probably  a  result  of  the  02^  A)  becoming  the 
dominant  species  in  the  afterglow  (rather  than  a  decrease  in  the  ozone 
production).  This  increased  ozone  deactivation  rate  in  the  afterglow 
causes  the  decrease  in  the  ozone  yields  at  the  detector  set  up. 

8.  Conclusions 

In  conclusion,  we  have  found  that  direct  nuclear  pumping  favors  the 
production  if  ozone.  In  addition  the  (^(^a)  yield  scales  with  the 
noble  gas  pressure  but  the  ozone  yield  scales  much  faster  for  a  given  power 
input.  Thus,  we  have  identified  the  photolysis  of  ozone  into  O2 ( ^A) 
as  a  very  promising  approach  for  the  efficient  production  of  large 
(^(^A)  yields.  (^(^A)  pressures  of  up  to  10  Torr  could  be  generated 
and  the  ozone  yields  necessary  for  this  scheme  have  indeed  been  observed 
experimentally. 

The  nuclear  sustained  electrical  discharge  generated  O2 ( ^A)  yields 
which  scaled  according  to  the  square  root  of  the  power.  No  saturation  was 
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28  Decrease  In  03/02(  A)  ratio  observed  as  a  function 
of  input  current.  This  graph  demonstrates  how  the 
low  E/N  discharge  favors  the  production  of  OgPA). 
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observed  up  to  the  limit  of  the  power  supply  used  which  resulted  in  an 
02(1  A)  concentration  of  1.4  mTorr  in  1  Torr  of  03-  In  addition, 
despite  an  increase  in  power  input  only  02(*a)  was  observed  to  increase 
substantially,  O3  yields  varied  only  slightly.  Finally,  the  diagnostics 
developed  for  this  work  made  the  very  sensitive  concentration  measurements 
possible  despite  the  adverse  conditions  caused  by  extraneous  radiation  from 
the  reactor. 
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KINETIC  MODEL  OF  PHOTOLYTIC  02(]A)  GENERATION 

ABSTRACT:  Kinetic  analysis  of  photolysis  of  03:02:He  mixtures  by  UV  light 
has  been  studied.  Dependence  of  02 ( ^ A )  yield  on  partial  pressures,  flash- 
lamp  intensity  and  duration  was  investigated.  Results  indicate  sufficient 
02(^A)  yield  for  operation  of  a  high-power  atomic  iodine  laser  at  the 
I*(2Pl/2)  -  ^-i(2p  )  transition. 

I.  Introduction 

2 

Operation  of  an  atomic  iodine  transfer  laser  at  the  I*(  Py2) 

1 -31yjq^|  ^2  ^  spin-orbit  transition  is  well  demonstrated  and  represents 

an  important  high-power  laser,  (Refs.  A1-A3) .  In  this  laser,  excited 

molecular  oxygen,  02( 1  A) ,  is  generated  chemically  and  then  transferred  to 

a  laser  cell  where  it  is  mixed  with  I This  results  in  the  dissociation 

2 

of  I2  and  the  subsequent  formation  of  '*<' Pl/2>  by  the  fast  near  resonant 
energy  transfer  reaction: 

o2(]A)  +  i(2p3/2)^i=±o2(3e)  +  i*(2p1/2).  (A1> 

The  equilibrium  constant  for  this  reaction  is  2.9  at  300°K.  Therefore,  in 

order  to  achieve  a  population  inversion,  it  is  necessary  that  the  ^2^  ratio 

[02(3Z)j 

be  greater  than  .17  at  room  temperature  (Ref.  A2). 


The  idea  of  achieving  a  population  inversion  on  th'e  atomic  iodine 

I(  p-j )  -*■  I(  P3/2)  Transition  by  near  resonant  energy  transfer  from 

O^A)  was  proposed  several  years  ago  (Ref.  A4).  The  lack  of  an  02(1A) 

generator  capable  of  delivering  a  sufficient  yield  prevented  the  rapid 
development  of  this  idea.  The  first  successful  operation  of  an  02(1a)-I2 

transfer  laser  used  a  chemical  source  for  the  generation  of  02(1a)  (Ref.  Al). 

In  this  laser  02(^A)  was  produced  by  reaction  of  chlorine  with  basic 
hydrogen  peroxide.  About  1  Torr  02  was  produced  with  an  02(V)  yield 
greater  than  40%. 

Although  operation  of  the  chemical  02(1a)-I2  laser  is  well  demonstrated 

and  has  proven  high  power  capabilities,  the  chemical  generation  of  02(1a) 

suffers  from  a  great  deal  of  shortcomings.  This  includes  the  inability  to 
generate  high  O^A)  density,  the  need  for  cooling,  and  the  use  of  corro¬ 
sive  materials.  Investigation  of  other  possible  sources  of  02 ( ^ A )  includes 
electron  beam  sustained  discharges  (Ref.  A5),  microwave  discharges  (Ref. 

A6),  direct  nuclear  pumping  (Ref.  A7)  and  photolytic  decomposition  of 
ozone  (Ref.  A8).  Recently  the  authors  described  the  prospects  of  utilizing 
nuclear  energy  to  simultaneously  produce  ozone  while  providing  energy  to  a 
KrF  nuclear  pumped  flashlamp  (X  =  249  nm)  for  decomposing  the  nuclear- 
produced  ozone  into  02(1A)  (Ref.  A9).  This  paper  describes  the  physical 
process  involved  in  02(*A)  production  by  photolysis  of  0g:02:He  mixtures. 

II.  Kinetic  Model 

In  the  photolysis  of  ozone  (220-300  nm  Hartley  band)  the  main 
dissociation  channels  are  (Ref.  A10): 
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03  +  hv  (220  <  X  <  300nm)  -  02(1A)  +  0(]D)  (A2) 

-  02(3Z)  +  0(3p) .  (A3) 

The  quantum  efficiency  for  channel  A2  is  determined  to  be  approximately 
90%  with  the  relative  quantum  yield  for  reactions  A2  and  A3  being 
constant  in  the  range  250  to  300  nm  (Ref.  A10).  For  X  >  310  nm  the 
primary  photodissociation  channel  is  (Ref.  All): 

03  +  hv  (X  >  300nm)  -*•  02(3Z,  ]A,  or  ]Z)  +  0(3p).  (A4) 

However,  the  photoabsorption  cross  section  of  ozone  in  this  range  is  much 
lower  than  for  X  <  300nm  (Fig.  1). 

The  rate  of  photodissociation  of  ozone  by  UV  light  extending  from 
X-j  to  X2  is  given  by 

H  • 

W  =  C  j  n(X)  a  (X)  dX  (A5) 

X1 

where  C  is  the  speed  of  light 

n(X)  is  the  density  of  photons  with  wavelength  X  ,  and 
a(X)  is  the  ozone  photoabsorption  cross  section  at  X. 

For  narrow-band  UV  light,  as  in  the  case  of  a  rare  gas  halide  excimer 

O 

flashlamp  (typical  bandwidth  =  20A),  a  good  approximation  of  W  is  given  by 

W  =  C  o  (X  )n 
o 


(A6) 


where  a(AQ)  is  the  absorption  cross  section  of  ozone  at  the  central 

wavelength  of  the  excimer  fluorescence  band, 
and  n  is  the  photon  density  in  the  fluorescence  band. 

Some  characteristics  of  several  excimer  flashlamps  that  can  be  utilized 
for  o  Ol)  generation  are  listed  in  Table  I. 

Table  A1 

Some  Properties  of  Excimer  Flashlamps  (Refs.  A12,  A13) 


System 

Wavelength, 

Ac  [nm] 

Maximum  Radiative 
Efficiency  [%] 

Ozone  Absorption  Cross 
Section  at  A0,  a[Mb] 

Xel* 

253 

37.2 

11.4 

KrF* 

249 

34.0 

11.1 

XeBr* 

282 

29.1 

3.4 

ArF* 

284 

22.4 

2.4 

The  use  of  a  standard  xenon  flashlamp  for  the  photolysis  of  ozone  is 
also  possible  at  a  reduced  efficiency  since  the  flashlamp  output  in  the 
220-300  nm  band  is  typically  less  than  10*  for  a  1  ms  pulse  and  decreases 
with  an  increase  in  the  flashlamp  duration  (Ref.  A14).  Although  UV  light 
above  310  nm  will  dissociate  ozone  mainly  into  02(  Z),  the  low 
photoabsorption  cross  section  of  ozone  in  this  region  limits  this  process. 

In  this  analysis  a  KrF*  flashlamp  is  considered  as  the  UV  source. 

Such  a  flashlamp  has  a  proven  efficiency  higher  than  10%  for  an  electrical 
discharge  and  electron  beam  pumping  (Ref.  A15),  and  an  estimated 


efficiency  of  14.25%  for  direct  nuclear  pumping  (Ref.  A13).  If  all  of  the 
flashlamp  output  is  at  the  central  wavelength  (X  =  249  nm),  then  the  photon 

density  in  the  absorption  cell  is  given  by: 

n  =  3.79  x  107I  cm-3  (A7) 

where  I  is  the  output  flashlamp  intensity  in  W/cm^. 

Therefore,  the  rate  of  ozone  dissociation  is  given  by 

W  =  13.88  I  sec'1.  (A8) 

Since  the  quantum  yield  of  reaction  (la)  is  90%,  the  rate  of  ozone  disso¬ 
ciation  through  channels  (la)  and  (lb)  is  given  by: 


Wla  = 

12.49  I 

sec”1 

(A9) 

“lb  = 

1.39  I 

sec'1 

(A10) 

The  reactions  considered  in  the  model  are  listed  in  Table  II.  Most  of  the 
rate  coefficients  considered  in  this  model  are  the  preferred  values  sug¬ 
gested  in  Reference  A12.  Generation  of  Og^A)  by  collisional  deactivation 
of  I)  by  ozone  was  neglected  since  the  branching  ratio  between 
reactions  (7b)  and  (7c)  are  unknown.  As  the  concentration  of  Ogt1!)  is  low, 
inclusion  of  reaction  (7c)  will  only  cause  a  slight  increase  in  the  produced 
O2 ( ) .  Reaction  17b  was  determined  to  proceed  £  3  times  faster  than 
reaction  17c;  therefore,  energy  transfer  from  O^D)  to  ground  state  molecular 
oxygen,  02(  I),  generating  02(  A)  was  neglected  in  this  model. 
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Table  A2 


Photolytic  Generation  of  02(^A)  Model 


Reaction 

No. 

Reaction 

Rate 

Coefficient3 

Refer¬ 

ence 

Photodecomposition  of  ozone: 

la 

03  +  hv  -  02(]A)  +  0^0) 

W-j  (see  text) 

lb 

-  02(3i)  +  0(3p) 

(see  text) 

Regeneration  of  ozone: 

2 

0(3p)  +  o2(3e)  +  o2(3z)  -  o3  +  o2(3z) 

6.95  x  10'34 

A12 

3 

0(3p)  +  02(3E)  +  He  03  +  He 

3.48  x  10"34 

A 16 

3 

Three-body  recombination  forming  02(  Z): 

4 

0(3p)  +  0(3p)  +  o2(3e)  ->  o2(3e)  +  o2(3e) 

7.2  x  10"33 

A17 

5 

0(3p)  +  0(3P)  +  He  ->  02(3Z)  +  He 

9.0  x  10'34 

A18 

Dissociation  of  ozone: 

6a 

03  +  0 ( 1 D )  -  02(3Z)  +  0(3P)  +  0(3p) 

1.2  x  10’10 

A 12 

6b 

-  o2(3e)  +  o2(3e) 

1.2  x  10"10 

A 12 

7a 

o3  +  O^Z)-  o2(3e)  +  o2(3z)  +  0(3p) 

1.54  x  10" 11 

A 12 

8 

V 

03  +  0(3p)  -  02(3E)  +  02(3E) 

-15 

8.4  x  10  13 

A 12 

9 

03  +  02(]a)  -  02(3e)  +  o2(3z)  +  0(3p) 

-15 

3.8  x  10 

A12 

aaflcaai^^ 

a  .sVijA  .  j  -.i 

' '  y\' 

Table  II  (Continued) 


Reaction 

No. 

Reaction 

Rate 

Coefficient3 

Refer¬ 

ence 

Reactions  deactivating  02(^Z)  in  addition 

to  7a: 

7b 

°3  +  02(1Z)  -  02(3I)  +  03 

-1 2 

6.6  x  10  xc 

A12 

7c 

-  o2(1a)  +  o3 

neglected 

A 12 

10 

o2(1i) 

+  0(3p)  -  o2(3z)  +  0(3p) 

8.0  x  10'14 

A12 

11 

o2(1z) 

+  o2(3i)  +  o2(3z)  +  o2(3r) 

4.0  x  10"17 

A 12 

12 

o2(1z) 

+  He  -  02(32)  +  He 

1.0  x  10'17 

A 12 

Reactions  deactivating  02(^A)  in  addition 

to  9: 

13 

C^A) 

+  0(3p)  -  o2(3e)  +  0(3p) 

1.3  x  10"16 

A19 

14 

o2(1a) 

+  02(1a)  -  Og^z)  +  02(3Z) 

-17 

2.0  x  10  17 

A20 

15 

O^A) 

+  02(3Z)  02(3E)  +  02(3Z) 

2.2  x  10-18 

A12 

16 

O^A) 

+  He  -  02(3E)  +  He 

1.0  x  10‘20 

A21 

Reactions  deactivating  0(^D)  in  addition 

to  6(a,b): 

17a 

0 ( 1 D)  • 

►  02(3Z)  -  0(3p)  +  02(]z) 

3.2  x  10"11 

A12 

17b 

-  0(3p  +  02(3z) 

-12 

8.0  x  10 

A12 

17c 

-  0(3p)  +  Q2(\) 

neglected 

A 12 

18 

0(1D)  • 

)■  He  -  0(3p)  +  He 

-15 

1.5  x  10  3 

A22 

aAll  rate  coefficients  are  in  units  of  molecules  cnr  at  298  K. 


*  -  •  V 


For  a  homogeneous  gas  mixture  consisting  of  o^one,  oxygen,  and  helium, 
and  assuming  the  mixture  remains  near  room  temperature  during  the  photolysis 


process,  the  rate  equations  governing  the  concentrations  of  the  different 
species  can  be  written  in  the  form: 


dN 


x  1  j 

-  I  (rate  of  reactions  creating  Nv)  -  t  (rate  of  reacti 


ons 


consuming  Nx) 


(All) 


where  the  rates  of  two-body  reactions  are  given  by  K12  [N^H^]  , 

the  rates  of  three-body  reactions  are  given  by  [N^ ][N2][N3] . 
and  the  rates  of  ozone  photodecomposition  are  given  by  w-j a [O3D 
for  reactions  la  and  Wi^tO^]  for  reaction  lb. 


Note  that  Wlfl  and  Wlb  depend  on  the  flash! amp  intensity  and  therefore 
follow  the  time  evolution  of  the  flashlamp  pulse.  In  this  model  W^a  and 
W^b  are  assumed  constant  throughout  the  absorption  cell  at  any  specific  time 
Physically,  this  means  that  the  photon  density  is  constant  in  the  absorp¬ 
tion  cell.  This  condition  is  valid  if  aCO^]  L  «  1  according  to  the  Beers- 

-18  2 

Lambert  law,  where  o  *  11.1  (10  )  cm  at  X  =  249  nm  and  L  is  the  width 

of  the  absorption  cell. 

For  example,  with  2  Torr  O3  this  condition  requires  L  «  1.33  cm  if 
the  photodecomposition  of  ozone  is  not  accounted  for.  When  the  cell  is 
irradiated  from  two  sides,  L  <<  2.66  %  .5  cm  would  satisfy  the  assumption 
of  uniform  photon  density  in  the  absorption  cell.  For  a  wider  cell  con¬ 
taining  >  2  Torr  0^,  it  would  be  necessary  to  include  the  effect  of  the 
spatial  photon  flux  reduction  as  it  propagates  inside  the  cell. 
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III.  Results  of  Calculations  and  Discussions 

A  typical  output  of  this  model  is  shown  in  Fig.  A2  which  demonstrates 

the  basic  kinetic  processes.  The  reactions  are  initiated  by  the  photo- 

11  3 

decomposition  of  O3  to  02(  A)  and  0(  D)  (reaction  la),  or  to  02(  Z)  and 

0(  p)  (reaction  lb).  The  branching  ratio  between  la  and  lb  is  consequently 

assumed  to  be  9:1.  The  flashlamp  intensity  as  a  function  of  time  is 

assumed  to  be  a  Gaussian  peaking  at  time  =  .042  seconds  with  a  period  of 

13.8  ms.  The  formed  0(^D)  is  quickly  deactivated  to  ground  state  atomic 

3  3  3 

oxygen,  0(  p),  by  collisions  with  03,  02  (  Z),  and  He.  The  0(  p)  then  re- 

combines  with  ground  state  molecular  oxygen,  02(  Z),  in  a  three-body  recom 

bi nation  reaction  to  form  ozone  (reactions  2  and  3).  The  reformation  of 

ozone  throughout  the  pulse  is  enhanced  at  high  He  pressures  and  is  responsible 

for  the  large  02( ^ A)  yields.  In  addition,  reactions  2  and  3  aid  in  the 

3 

depletion  of  the  02(  Z)  population  which  further  enhances  the 

[02(’a)]  [o,('a)] 

— =-r= - ratio.  Large - 5 - ratios  make  the  forward  path  of  reaction  A1 

[o2(3z)]  [o2(3z)] 

2 

more  favorable  and  thus  more  I*(  py2)  can  be  Produced.  As  time  progresses, 

the  02(1A)  concentration  is  damped  due  to  collisional  processes  while  the 

02(  Z)  density  exceeds  its  initial  value  as  a  result  of  the  decomposition 
of  the  ozone. 

Fig.  A3  shows  the  effect  of  helium  pressure  on  02(^A)  production.  The 
02(^A)  yield  increases  with  helium  pressure.  This  is  due  to  the  reformation 
of  ozone  throughout  the  pulse  by  three-body  recombination  (reaction  3), 
which  is  favored  at  higher  pressures.  Moreover,  since  the  rate  coefficient 


MAX  02(’A)  PARTIAL  PRESSURE 


aw 


FSL-83-17! 


Fig.  A3.  Maximum  02(1A)  partial  pressure  (T  *  298  K)  and 

Cop(1a)3 

- * —  ratio  as  a  function  of  He  pressure  for 

C02(3Z)] 

different  peak  flashlamp  intensities 
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for  the  deactivation  of  02(  A)  by  helium  is  low  (<10  ),  such  high 

helium  pressures  do  not  greatly  reduce  the  effective  lifetime  of  Og ( ^ A ) . 

The  three-body  recombination  reforming  02(  E)  (reaction  5)  competes  with 

this  positive  feedback  regeneration  of  03  and  limits  its  effect.  As  the 

helium  pressure  is  further  increased,  the  concentration  of  02(1A) 
reaches  a  saturation  point.  This  saturation  point  is  due  to  the  quench¬ 
ing  of  O-^A)  by  the  regenerated  03  (reaction  6)  which  competes  with  its 
photodecomposition  into  02(^A).  The  saturation  point  of  02 ( ^ A )  yield  occurs 
at  lower  helium  pressures  as  the  flashlamp  intensity  is  decreased.  This 
is  an  expected  effect  of  the  reduction  in  the  rate  of  0^  photodecomposition. 

The  production  efficiency  of  02(^A)  is  highest  for  a  short  flashlamp 
duration.  Fig.  A4is  a  plot  of  the  maximum  02(1A)  yield  as  a  function  of 
pulse  width  (e-^  of  peak  value)  for  a  constant  peak  flashlamp  intensity  and 
different  02:03  mixtures.  As  the  pulse  duration  is  increased  the  02(^A) 

yield  increases  up  to  the  steady  state  solution  of  this  model.  However, 
the  02(  A)  yield  for  a  given  absorbed  energy  is  reduced  with  the  longer 
pulses.  This  reduction  in  02(^A)  production  efficiency  is  greater  for 
large  0^  fractions  because  it  is  a  result  of  the  slow  decomposition  of  03 

which  quenches  02(^A). 

For  operation  of  an  02(^A)-I2  laser  it  is  necessary  to  transfer  02(^A) 

from  its  generation  point  to  the  laser  cell;  therefore,  the  laser  effi¬ 
ciency  will  depend  on  the  effective  lifetime  of  02(^A).  Fig.  A5  shows  the 
time  behavior  of  02(^A)  and  0^  concentrations  for  different  peak  flashlamp 


intensities.  The  02(^A)  yield  and  its  effective  lifetime  increase  with 
increasing  flashlamp  intensity  due  to  the  enhanced  rate  of  photodecomposi¬ 
tion  of  O^.  For  2  Torr  0^.  18  Torr  02,  3800  Torr  He,  a  peak  flashlamp 

2 

intensity  of  1  KW/cm  ,  and  a  pulse  width  of  13.8  ms,  the  maximum  partial 

O2 ( ^ A )  pressure  is  7.5  Torr  with  an  j ]  ratio  of  .58  and  an  effective  02(^A) 

halflife  of  ^  70  ms.  Further  enhancement  of  the  effective  halflife  of 
02(^A)  is  possible  using  a  more  intense  flashlamp  to  further  decompose  O^. 


IV.  Conclusion 


Kinetic  analysis  of  photolysis  of  several  atmospheres  of  He 

mixtures  predicts  sufficient  02(^A)  yield  for  the  operation  of  an 

02(^A)-I2  transfer  laser.  The  enhancement  of  the  rate  of  three-body 
regeneration  of  ozone  under  high  pressures  is  responsible  for  the  large 
02(1A)  yield.  Applications  of  such  a  system  can  include  nuclear  pumping  of 
an  iodine  laser  by  utilizing  a  nucl ear-pumped  fluorescence  source  or  its 
development  as  a  solar-pumped  transfer  laser.  In  addition,  the  photolytic 
pumped  02(^A)-I2  laser  offers  the  potential  of  a  continuous  high-power  laser 
with  complete  chemical  recirculation. 
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Appendix  B. 

DISCUSSION  OF  THE  MODEL 

1)  Nuclear  Pumping 

The  basic  processes  modeled  include  the  ionization  and  excitation 
of  argon  atoms  and  the  oxygen  molecules.  The  latter  has  been 
calculated  with  an  electron  distribution  function  and  will  be 
discussed  in  a  later  section.  The  ionization  of  argon  and  oxygen  can 
be  calculated  using  the  W  value  approach,  where  the  W  value  is  the 
amount  of  energy  required  to  create  a  single  ion  pair.  Therefore,  the 
source  rate  of  ions  becomes  (Eq.  B3) : 

where,  Edep  =  normalized  energy  deposition  (figure  Bl) , 

4>n  =  thermal  neutron  flux, 

W  =  26.4  ev/i-p  for  Ar,  and  (Ref.  B2) 

32.5  ev/i-p  for  02. 

This  is  also  used  as  the  source  rate  for  the  fast  electrons  in  the 
model . 

A  similar  approach  has  been  used  to  describe  the  source  rate  of 
Ar  metastables  with  an  effective  W  value  (Eq.  B4)  (Ref.  B3) : 

(4)  S*  =  S+ 

which  gave  an  Ar  concentrations  similar  to  the  electron  distribution 
function  approach  (Eq.  B5)  (Ref.  B4) : 

(5)  Ar  ♦  e"  t  Ar*  +  e"  Kf  =  3 . 1  ( 1 0” 1 3 )  cm3/sec 

Kr  =  1 .7(1 0~ 9 )  cm3/sec 

Therefore,  the  nuclear  pumping  portion  of  the  model  consists  of 
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Figure  B1  Power  deposition  of  Ions  vs.  radial  distance  for  the  case  of  helium 
gas  with  boron  coatlnq. 

A  —  Chung's  result,  pressure  at  760  torr 

B  —  Guyot's  result,  pressure  at  760  torr 

C  —  Chung's  result,  pressure  at  100  torr 

D  — -  Guyot's  result,  pressure  at  100  torr 

Coating  thickness  is  1.739  um 


using  the  appropriate  power  deposition  term  and  coupling  the  neutron 
flux  to  the  model  with  equation  B3.  Several  power  deposition 
calculations  have  been  performed  and  the  one  quoted  here  is  from  the 
most  recent  study  (Refs.  B5,  B6) .  The  effects  of  the  spatial 
distribution  of  the  power  (figure  Bl)  are  discussed  later. 

2)  Argon  Model 

Nuclear  pumping  of  an  Ar-02  mixture  can  be  broken  down  into;  1) 
Nuclear  pumping  of  Ar,  2)  energy  transfer  from  Ar  to  0 2,  and  3)  energy 
redistribution  in  0^.  The  net  effect  of  adding  a  diatomic  specie,  such 
as  oxygen,  to  a  noble  gas  discharge  is  the  complete  transfer  of  the 
energy  deposited  in  it  to  the  diatomic. 

The  reactions  of  the  argon  system  can  be  found  in  Table  Bl.  A 
flow  diagram  for  these  reactions  is  shown  in  figure  B2  and  reveals  a 
very  efficient  pathway  for  the  formation  of  0*  ( 96%  of  the  energy 

deposited  in  Ar+  and  Ar*).  This  efficient  energy  transfer  is 

responsible  for  the  quenching  of  the  optical  emissions  from  Ar2  and 

* 

Ar  .  The  charge  transfer  reaction  (Eq.  B6) : 

(6)  Ar*  +  02  02  +  Ar  +  Ar, 

prevents  the  dissociative  recombination  reaction  (Eq.  B7)  : 

(7)  Ar*  +  e~  Ar  +  Ar  , 

which  effectively  decreases  the  Ar*  densities.  Therefore,  if  the 
Ar  densities  are  suppressed  then,  the  reactions  (Eqs.  B8,  B9) : 

(8)  Ar  +  Ar  -*■  Ar2, 

(9)  Ar2  -+■  Ar  +  Ar  +  hv, 

are  prohibited.  Consequently,  almost  all  optical  emissions  from  these 
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Table  B1 


ARGON  MODEL  REACTIONS  (Eqs.  B10-B21) 
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i»tl  ~  ■  — -m  ~  w  -  w  r*  <*  r  *  "•»  —  m  r- *»  -•>  ' 
i  *  «  «■  «  4  *  ,  r,  •  ,  *  . 


cm 


/sec 


Reference 
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Ar  +  a  -*■  Ar+  +  e-  +  a 

K  = 

1.3 
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B2,  B5 
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* 
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4.1 
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B3 
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Figure  B2  Energy  flow  diagram  in  an  Argon-0 
nuclear  pumped  plasma. 

Net  efficiency:  a  -*■  Ar  09:  96 


states  will  be  suppressed  upon  adding  even  an  impurity  of  0 ^  to  the 
system . 

3)  Oxygen  Model 

The  reaction  set  listed  in  Table  B2  is  based  on  a  two  energy  group 
electron  distribution  function.  The  energy  groups  were  defined  as: 


for 

10“2  <  E 

< 

leV 

for 

E 

> 

leV 

This  division  proved  to  be  convenient  when  classifying  the  vacious 
superelastic  collisions  into  the  cases  which  could  supply  sufficient 
energy  to  the  electrons  to  move  them  from  group  2  to  group  1  and  those 
that  couldn't.  In  addition,  the  four  distribution  functions  shown  in 
figure  3  were  observed  by  each  author  to  have  a  low  energy  portion 

(T  a  .1  eV)  which  was  independent  of  the  high  energy  tail. 

6 

A  comparison  of  the  electron  energy  degradation  spectrum  for  four 
cases  is  shown  in  figure  B3.  The  results  of  Lo  and  Miley  (Ref.  BIO), 
have  often  been  accepted  as  describing  the  energy  distribution  function 
for  1  MeV  alpha  particles  in  helium.  However,  the  addition  of  oxygen 
to  the  mixture  should  cause  a  downshift  in  the  spectrum.  Therefore, 
the  use  of  this  spectrum  to  calculate  excitation  rates  from  the  cross 
sections  of  Fournier  (Ref.  Bll)  (figure  B4)  is  probably  an  overestima¬ 
tion  (Table  B2).  A  better  approximation  is  the  result  of  Fournier 
(Ref.  Bll)  for  an  Ar-t^  mixture  at  1  atmosphere  and  pumped  by  a  2  KeV 
e-beam.  The  reaction  rates  calculated  using  this  result  for  E  >  1 
eV  are  listed  in  table  B2  along  with  the  He  case.  The  technique  used 
by  Fournier  to  calculate  the  distribution  functions  was  benchmarked  by 
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Table  B2 

REACTION  RATES  (Eqs  .  B22-B66) 


<  <7v> 


Primary  Processes  (excl.atlon) 

[cm  /sec) 
Ct-partlcle 

lcm3secl 

2  keV  e-beam 

Eqn.  No, 

(He) 

(Ar-OJ 

2 

«2  ♦  °2  <x^  >  -►  e~  +  °2(a  'A) 

1.05  (-15) 

1.05  (-15) 

B  22 

e”  +  02  (X^)  f  +  02<a'A) 

3.74  (-12) 

9.65  (-13) 

B23 

e’  ♦  0,  (X^)  -*•  e’  +  0  Vi) 

2  2  2  2 

1.75  (-19) 

1.75  (-19) 

B2  4 

e~  +  02  (X3Z>  -*■  e”  +  °2tt>,^) 

9.99  (-13) 

2.86  (-13) 

B25 

e”  +  02  (X^!)  -*■  e2  +  02  <A3D 

3.60  (-12) 

1.75  (-13) 

B2  6 

e~  ♦  02  <X3I)  «2  +  °2  <b3^ 

6.37  (-13) 

7.75  (-14) 

B27 

*2  +  °2  <a'^  ■“  *2  +  °2 

1.71  (-12) 

1.71  (-12) 

B28 

e”  +  02  (a'^J  "*■  *2  +  °2 

2.57  (-11) 

8.56  (-12) 

B2  9 

*1  *  °2  <a<^>  *2  *  °2  <a3“J 

3.60  (-12) 

1.75  (-13) 

B3  0 

♦  02  (a'fi)  "*■  *2  °2  (b3^') 

6.37  (-13) 

7.75  (-14) 

B3  1 

e~  ♦  02  (b’E)  *2  +  °2  <a3i) 

3.6  (-12) 

1.75  (-13) 

B3  2 

*1  +  °2  *2  +  °2  (b3^J 

6.37  (-13) 

7.75  (-14) 

B33 

e~  ♦  02  (A3!)  -v  e”  ♦  02  (B3!) 

8.1  (-12) 

8.1  (-13) 

B34 

Table  B2  (Continued) 

REACTION  RATES  (Eqs.  B22-B66) 
<  o  v> 


m 

ll 


(cm  /seel 

Primary  Processes  (Superelastic)  Cl-particle 


2  '  2  ’1  ~2  '  “ 

02  (B3S)  +  e^  e“  +  02  (X3!) 

02  <A3D  +  ®2*K  +  °2 

02  (A3E)  ♦  e”  -*■  e ~  ♦  02  (X3E) 

0„  (b'E)  +  e"  e”  +  0,  (X3£) 

2  2  2  2 

02  (b’E)  +  e~  •*  e"  +  02  (X3£) 

0„  (a’A)  +  e“-*-  e”  +  0,  (X^) 
2  2  2  2 


02  (a'A)  ♦  9~  •*  9~  *  02  (X%) 


2  t  2 


1  1  2 


02  (a!  )  +  *2  +  °2  <8*A> 


1  1  2 


2  2  2 


1  1  2 


2  1  2 


I  1  2 


2  1  2 


1  ’1  2 


2  2  2  2 
02  (B3^)  ♦  ej  -*•  e”  +  02 


2.76 

(-6) 

1.96 

(-10) 

6.46 

(-9) 

2.73 

(-10) 

3.02 

(-10) 

7.67 

(-12) 

9.63 

(-11) 

1.08 

(-11) 

1.84 

(-8) 

1.30 

(-10) 

8.48 

(-9) 

1.82 

(-10) 

2.51 

(-9) 

6.63 

(-11) 

9.20 

(-9) 

6.50 

(-11) 

2.83 

(-9) 

9.1  ( 

-11) 

1.22 

(-8) 

1.08 

(-10) 

y.V.V  ‘i  -  aV-'-'-V.  i*-  «*-  <*i  i' 


1cm3 seel 

2  keV  e-beam 

m 

(Ar-02) 

2.76  (-6) 

B3  5 

7.25  (-11) 

B36 

8.48  (-9) 

B37 

8.26  (-11) 

B3B 

3.02  (-10) 

B3  9 

2.77  (-12) 

B4  0 

9.83  (-11) 

B4  1 

3.42  (-12) 

B4  2 

1.84  (-8) 

B4  3 

4.83  (-11) 

B44 

8.48  (-9) 

B4  5 

5.51  (-11) 

B46 

2.51  (-9) 

B  4  7 

2.34  (-11) 

B48 

4.37  (-9) 

B49 

2.4  (-11) 

B  5  0 

1.70  (-9) 

B  5  1 

2.75  (-11) 

B52 

1.22  (-8) 

B  5  3 

3.5  (-11) 

B  5  4 

•  •  ■  •  ’  .  '  V  v'  •  -  •  ■ 

Table  B2  (Concluded) 

REACTION  rates  (Eqs  .  B22-B66) 

<  <Jv> _  _ 


I cm  /seel 

lcm3secl 

■gl 

Primary  Processes  (Ionization) 

Ot-oartlcle 

2  keV  e-beam 

(He) 

(Ar-02) 

02  (X%>.  (aIA),(b,E),(A32),(B32:) 

+  e,  +  e  +  0* 

1  2  2 

8.68  (-13) 

1.48  (-13) 

B56 

—  + 

-^•2  ♦  0  *4“  0  +  0 

1.26  (-13) 

3.40  (-14) 

B  5  6 

Dissociative  Attachment 

02  (X3!)  +  e”-*-  0  +  o” 

1.20  (-13) 

9.96  (-15) 

B57 

0  (a%  +  e'-*-  0  +  0~ 

6.6  (-13) 

8.7  (-14) 

B5  8 

02  (b'l)  +  e^  0  +  o" 

9.05  (-13) 

1.28  (-13) 

B  5  9 

02  (A3!)  +  e~  -*  0  ♦  0~ 

9.5  (-13) 

3.15  (-13) 

B6  0 

0,  (A^E)  +  «“-*•  0  ♦  0* 

2  2 

1.05  (-15) 

1.05  (-15) 

B6  1 

02  (B3E)  +  •' -t  0  +  0~ 

1.24  (-12) 

4.58  (-13) 

B6  2 

0,  (B3E)  +  e"  -►  0  +  o" 

2  2 

2.31  (-14) 

2.31  (-14) 

B63 

4.5  eV  Loss  Processes 

0,  (X3!)  ♦  e’  -»-e’  +  0* 

2  12  2 

8.02  (-13) 

5.83  (-14) 

B6  4 

Elastic  Col  1 Isions 

02  (X%)  +  e"-t  a”  ♦  02 

8.00  (-10) 

2.84  (-10) 

B65 

02  (X3E)  +  e"-*  e'  +  °2 

2.14  (-8) 

B6  6 
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performing  a  calculation  in  ^  and  comparing  it  to  the  distribution 
function  measured  by  Shure  (Ref.  B12)  in  a  hollow  cathode  discharge. 

In  general,  the  agreement  was  good  over  the  energy  range  examined 
(  1  <  E  <  10  eV) . 

The  calculation  of  <ov>  was  performed  with  cross  sections 
obtained  from  Fournier  (Ref.  Bll).  In  the  model,  the  following 
processes  have  been  considered  (table  B2) : 

(A)  Ionization; 

02(X3!:),  Og^A),  02(b1Z),  02(A3I),  and  02(B3E)  , 

(B)  Dissociative  ionization  of  the  same  states, 

(C)  Dissociative  attachment  of  the  same  states, 

(D)  Excitation  reactions; 

02(X3!)  ^02(a1A),02(b1E),02(A3E),  and  02(B3E), 

OgU’A)  -*  02(b1I),  02(A3E>,  and  02(B3E), 

O^b1!)  -  02(A3E)  and  02(B3E), 

02(A3E)  ■-  02(B3E), 

(E)  Superelastic  collisions; 

02(B3E)  ->02(A3E), Ogt^E)  ^(a^),  and  02(X3E), 

02(A3E)  02(b1E),  OgCa^),  and  02(X3E) , 

Og^E)-  Og^A)  and  02(X3E), 

OgCa’A)  ■*»  02(X3E). 

The  cross  section  for  many  of  these  processes  had  to  be  generated  from 
the  known  set  of  cross  sections  given  in  figure  B4.  To  obtain  an 
excitation  or  ionization  rate  from  a  given  electronic  state,  the  cross 
section  for  the  ground  state  process  was  translated  to  a  lower  energy 
by  the  energy  content  of  the  electronic  state  of  interest.  If 
possible,  the  resulting  cross  section  was  normalized  to  an 


79 


experimental  results  ((^(a^A)  -^(^(b^E)  )  (Ref.  B13).  The  superelastic 
rates  were  then  calculated  from  the  excitation  cross  sections  generated 
from  the  ground  state  processes. 

The  results  of  these  calculations  are  listed  in  Table  B2.  It  is 

apparent  from  this  table,  that  there  is  a  strong  coupling  between  the 

Og(a  A),  02(b  E)  ,and  Og(X  E)  states.  Superelastic  rates  from  the 

Og(A3E)  and  Og(B3E)  states  to  the  lower  states  are  large  and  suggest  a 

strong  relaxation  of  the  higher  lying  states  during  the  discharge. 

3  3 

However,  the  excitation  of  the  Og(A  E)  and  Og(B  E)  by  electrons  is 
weak  in  the  two  cases  studied  here.  Consequently,  any  contributions 
to  the  0  (b^Z)  and  Og^A)  states  will  be  small  since  the  pumping 
rates  into  the  higher  states  are  so  low.  A  similar  argument  can  be 
applied  to  the  Og(c3A)  and  OgCe^E)  states  which  lie  near  the  Og(A3E) 
state. 

4)  Chemical  Kinetics 

Two  species  were  measured  in  the  experiments:  Og^A)  and  0^. 
Therefore,  this  section  is  subdivided  into  two  topics.  The  first 
discussion  describes  the  reactions  which  generate  O^^A)  while  the 
second  section  discusses  the  reactions  important  in  the  formation  of 
ozone . 

a.  02(a  a) 

The  development  of  the  Ogte^A)  model  began  by  considering  the  six 
pathways  capable  of  generating  OgCaV)  (Eqs.  B67-B72): 

(67)  e"  02  -*•  OgU^)  +  6g , 

(68)  0*  ♦  e2  OgU^)  +  e", 
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(69)  0(3P)  +  0(3P)  +  M  Og^A)  +  M, 

(70)  03  +  o(3P)  ♦  (^(a’A)  +  02, 

(71)  Og  +  Og  "*■  02(a^A)  ■*■  Og,  and 

(72) 0(1D)+02  -  OgU1*)  +  0(3P). 

Among  these  reactions,  70  and  71  were  found  to  be  the  least  important 
in  comparison  to  reactions  67  and  68.  Attempts  to  classify  reaction 
69  was  difficult,  because  the  rate  has  only  been  calculated  or 
estimated  in  the  research  performed  to  date.  This  is  apparent  in 
Table  B3  which  is  a  summary  of  the  current  literature  on  atomic  oxygen 
recombination.  In  this  table,  reaction  rate  73  is  inconsistent  with 
other  estimations.  Since  it  was  obtained  through  a  calculation  while 
the  others  were  based  on  experimental  observation,  the  other  estimates 
are  probably  better.  Thus,  based  on  the  estimates  available  (table 
3),  there  should  be  no  significant  contribution  to  the  (^(a1^) 
concentration  from  atomic  recombination.  Consequently,  the  dominant 
formation  pathways  in  a  nuclear  or  e-beam  pumped  system  appear  to  be 
reactions  67,  68  and  72. 

Of  the  five  states  above  Ogta^),  only  OgCb1!)  is  not  readily 
dissociated  by  electron  impact.  This  is  probably  due  to  its 
relatively  low  energy  compared  to  the  other  states  (1.62  eV  versus  E  > 
4.5  eV)  of  interest.  In  addition,  the  OgCb^E)  state  was  the  only  state 
which  was  able  to  build  up  a  significant  concentration  throughout  the 
pulse.  Consequently,  it  is  the  only  state  which  is  of  any  importance 
in  the  formation  of  02(a^A).  This  conclusion  is  supported  by  the 
reaction  rates  calculated  in  the  last  section  and  by  the  rates  listed 
in  Table  B3.  Even  though  the  recombination  rates  into  the  02(A3£)  and 
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(73) 

1.2  (10"JJ 

(74) 

7.4(10"33) 

(75) 

( 10"33) 

(76) 

2.8(10-33) 

(77) 

9.4(10"32) 

(78) 

- 

(79) 

(80) 

do"33) 

(81) 

“ 

(82) 

- 

(83) 

- 

(84) 

- 

(85) 

3(10“33) 

Not 

Favored 


Not 

Observable 


ado"37)1 


Table  B3  (Units  of  cm°/sec) (Eqs.  B73-B85) 
02(b1l)  I  0„(A3£)  I  Comments  I 


0„(A3r) 


3.2(10"34} 


3.8(10" 


Easily 

Dissociated 


1. 1( 10"37) 


1.7(10"37*1  3.3( 10"37) 


1 . 4( 10"37 ) 


Calc. 

Exp.  Disch. 
Calc. 

Exp. Shock  Tube 

Exp. Shock  Tube 

Exp.  No  +  N 
Exp.  No  +  N 

Exp.  No  +  N 
Exp.  Disch. 

Exp.  No  +  N 
Unknown 
Summary 


Two  body  dependence 
Estimated  from  02(^2)  behavior 


24 

25,  26 


V, 


82 


OgCB^Z)  states  are  the  largest,  the  majority  of  the  0  atoms  are  lost 
by  the  three  body  reaction  (Eq.  B86) : 

(86)  0  +  Og  +  Ar  ■+  O3  +  Ar, 

in  forming  0^.  This  effect  is  obvious  for  [02]>>[0]  since  reaction  69 
is  second  order  in  [0]  while  86  is  first  order. 

To  complete  the  model  of  02(a^A)  generation,  it  is  necessary  to 
examine  the  possible  02(b1Z)  production  mechanisms.  Again,  there  are 
several  candidates  (Eqs.  B87-B92): 

(87)  e"  +  02  02(b1I)  +  e", 

(88)  e2  +  02  •+  02(b^Z)  +  e2, 

(89)  e"  +  o2(z  )  -  02(b1Z)  +  e", 

(90)  0( ’d)  +  02  -»■  02(b1Z)  +  0(3P), 

(91)  0  +  0  +  Ar  +  02(b1Z)  +  Ar,  and 

(92)  0+  +  0"  -  02(b1Z)  +  02, 

where,  reaction  91  can  be  eliminated  with  the  arguments  presented 
above  and  reaction  92  was  found  to  be  negligible.  In  addition, 
reactions  87,88,  and  89  were  found  to  be  insignificant  in  comparison 
to  the  pumping  rate  of  reaction  90.  This  reaction  was  found  to 
dominated  both  the  generation  of  02(b^Z)  and  the  loss  of  O^D). 

The  main  source  of  O^D)  in  the  discharge  is  the  dissociative 
recombination  of  0*  (Eqs.  B93-B95): 

Branching 

Ratio 

(93)  02  +  e"  +0+0,  0.008 

(94)  -  0(1D)  +  0(1D),  0.025 

+  0(1S)  +  0(1S), 


(95) 


0.002 


.HUS 


V 
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where,  the  Og  is  also  lost  by  the  reactions  (Eqs.  B96-B99)  : 


(96)  0*  +  e2*  e2  °2  +  e2*  0.006 

(97)  Og  +  e^  ♦  Ar  Og  +  Ar,  0.013 

(98)  Og  +  02  +  Ar  oj  +  Ar,  0.93 

(99)  +  °g  oj  +  0g,  0.012 


and  the  0^  is  consumed  by  a  variety  of  fast  reactions  but  dominated  by 
dissociative  recombination.  The  final  model  for  the  production  of 
Og(a1A)  is  diagrammed  in  figure  B5»  In  this  figure,  the  relative 
importance  of  the  Ogfe^A)  pumping  mechanisms  are;  (67)  =  0.08, 
(68)  =  0.77,  and  (72)  =  0.14.  From  this  breakdown,  it  is  apparent 
that  the  main  excitation  mechanism  is  the  energy  transfer  from  0 ( 1D ) 
to  Og  to  form  Og(b1J:)  followed  by  the  relaxation  of  OgCb1^)  into 
Og^A). 

Throughout  the  discharge  and  afterglow,  several  of  the  species 
generated  can  lead  to  the  deactivation  of  Ogta1 A)  and  Og(b1Z).  A  study 
of  the  known  reactions  found  the  following  reactions  to  be  important 
losses  of  0  (a^A)  (Eqs.  B100-B107): 

Discharge 

(100)  Og(a^A)  +  0(1S)  -*•  0(  1D)  +  Ogtb1):), 

(101)  ♦  0+0+0, 

(102)  OgU’A)  +  0g  +  20g  +  e", 

(103)  Og^A)  +  eg  +  0g  +  eg, 

Afterglow 

(104)  Og(a^A)  +  Ar  •+  0g  +  Ar, 

(105)  Og^A)  +  Og  +  20g, 

(106)  Og^A)  +02-*-  20 g  +  0,  and 


Figure  B5  Energy, transfer  leading  to  the  formati 
of  Op (A)  dominates  all  of  the  other 
production  mechanisms  in  a  nuclear 
pumped  plasma. 


(107)  02(a'A)  +  wall  -*•  Og. 

At  the  low  OgCa^A)  densities  predicted,  the  dimolecular  reaction  was 
found  to  be  unimportant.  Similarly,  the  dominant  losses  for  O^b1!) 
are  (Eqs.  B108-B112): 

(108)  02(bh)  +  03  +  202  +  0, 

(109)  02(b^Z)  +  e2  -*■  Og^A)  +  e2> 

(110)  02(b1Z)  +  Ar  -*■  02  +  Ar, 

(111)  OgCb1!)  +  eg  02  +  e2,  and 

(112)  02(b1Z)  +  wall  -*■  02, 

where,  reaction  108  is  very  important  in  the  0^  kinetics  scheme. 


b.  Ozone 

During  this  study,  0^  was  found  to  be  one  of  the  major  species  in 
an  Ar-02  afterglow.  Both  calculations  and  experiments  have  indicated 
this  which  leads  to  the  present  discussion  of  the  kinetics  involved  in 
the  formation  of  0^.  As  was  the  case  in  the  last  section,  there  are 
many  potential  sources  of  0^  in  the  discharge.  With  the  aid  of  the 
model  it  was  possible  to  sort  through  each  of  the  0^  producing 
reactions  known  and  determine  their  relative  importance  (Eqs.  B113-B117) 

Contributing 

Fraction 


(113) 

0  +■  02  +  Ar  ■*  Ar, 

0.65 

(114) 

oj  +  o  o3  +  o2, 

0.22 

(115) 

0  +  O'  -  03  +  e”, 

0.08 

(116) 

0_  +  09  0,  +  0-, 

0.04 

© 


with  the  majority  of  the  0^  being  formed  by  the  three  body 
recombination  reaction  (113)*  The  source  of  0  atoms  and  0*  ean  be 
traced  to  02«  A  diagram  of  this  model  is  indicated  in  figure  B6  with 
the  branching  ratio  of  each  important  reaction. 

The  inherent  stability  of  0^  results  in  the  formation  of  large 
concentrations  (with  respect  to  OgCa1^))  after  the  discharge.  Indeed, 
only  two  loss  mechanisms  were  found  important  throughout  the  duration 
of  the  discharge  (Eqs.  B108-B109): 

(108)  03  +  q^S)  -  202, 

(109)  03  +  Ogtb1!)  ■>  202  +  0, 

but  since  reaction  (109)  results  in  the  formation  of  atomic  oxygen,  at 
least  1/3  of  the  0  atoms  will  return  to  form  03  under  the  conditions 
of  Ar  pressure  <  150  Torr.  If  the  Ar  pressure  »  150  torr,  then  the 
effects  of  diffusion  become  negligible  and  all  of  the  atomic  oxygen 
recombines  to  form  ozone.  This  will  result  in  reaction  (108)  becoming 
the  dominant  loss  of  ozone  throughout  the  pulse.  However, 
at  P(Ar)  =150  Torr,  the  relative  importance  of  these  two  reactions 
were:  (108)  =  0.45,  and  (109)  -  0.44.  Therefore,  any  loss  of  03  in  a 
long  pulse  can  be  attributed  to  the  reactions  which  occur  during  the 
collisions  with  oOs)  and  02(b1l).  Since,  these  states  are  a  result  of 
the  recombination  of  0*,  it  appears  that  a  short  pulse  will  be  more 
efficient  in  the  formation  of  03  because  the  metastable  species 
require  a  long  time  to  build  up. 
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One  of  the  largest  discrepancies  observed  so  far  between  the  modeling 
and  experiments  involves  ozone  yields.  The  production  of  ozone  is  thought 
to  be  due  to  the  three  body  reaction  (Eq.  B110): 

(110)  0  +  02  +  M  ■*>  03  +  M. 

Thus,  any  discrepancy  must  be  a  direct  result  of  an  overestimation  of  the 
atomic  oxygen  concentrations.  One  potential  reason  for  this  problem  may  be 
the  neglected  wall  interaction  mechanisms  in  the  model.  Attempts  to  model 
this  effect  immediately  produced  better  agreement  but  required  a  greatly 
enhanced  reaction  rate  in  comparison  to  the  value  obtained  by  the  general 
expression  (Eq.  Bill): 

(111)  Td  »  4  P  +  iy 

where: 

2  2 

A  ■  (R/tt)  *  characteristic  dimension 
P  *  pressure  of  gas  :  150  Torr 

2 

D  -  diffusion  coefficient:  228  cm  Torr/sec  for  0  atoms  (Ref.  B24) 

R  ■  radius  of  tube  ■  1.27  cm 

4 

G  ■  thermal  velocity  6.3  x  10  cm/sec 
y  ■  wall  deactivation  efficiency,  10  ^  <  y  <  1.0 
This  expression  yields  l/x^  ■  9  sec  ^  for  a  cosine  function  distribution. 
Consequently,  it  was  necessary  to  determine  more  details  about  the 
diffusion  mechanism  and  associated  spatial  effects. 

The  power  deposition  profile  is  similar  to  Figure  B2  which  is  a  result 
derived  by  A.  Chung  (Ref.  B5)  for  various  helium  pressures.  The  very  steep 
gradients  resolved  can  strongly  influence  diffusion.  To  obtain  an  estimate 
of  the  importance  of  this  profile  an  analytic  solution  was  obtained  to  the 
diffusion  equation.  The  source  was  represented  as  a  square  pulse  in  time 


with  the  following  spatial  distribution  (Eq.  B 1 12) : 

(112)  S  =  A  (1  -  sin  Bx) 

o  o 

where: 

17  3 

A  =  source  rate  of  0  atoms  ~  10  a/cm  -sec 
o 

B  =  tt/2. 

The  solution  to  the  diffusion  equation  (Eq.  B113): 

2 

(113)  =  D  d— A  (1  -  sin  Bx)  U(t-tQ) 

dt  dt 

subject  to  the  B.C.  (Eqs.  B114,  B 1 15 ) : 

(114)  «  o  @  X  -  1.0 
dx 

(115)  [0]  -  0  @  X  -  0 

was  easily  obtained  via  Laplace  transform  techniques.  Homogeneous  reactions 
were  then  added  to  the  solution  via  Danckwert's  method  (Eq.  B116) (Ref.  B27) . 

rc  "V  ,  ~V 

(116)  [0]  -  J  [0],  e  H  dt’  +  CO],  e 

o 

A  typical  result  for  a  10  millisecond  pulse  and  no  homogeneous  reactions 
is  shown  in  Figure  B7.  Similarly,  if  =  50  only  a  small  decrease  in  the 
peak  concentrations  is  observed  in  Figure  88  while  the  strong  depletion 
after  the  pulse  is  expected.  Consequently,  the  diffusion  on  this  simple 
slab  case  appears  to  dominate  the  losses  of  atomic  oxygen  during  the  pulse. 
This  is  further  exemplified  in  Table  84  which  lists  the  diffusion  decay 
constant  as  a  function  of  position  in  the  slab. 

A  comparison  of  the  decay  constants  @  t  =  .012  msec  and  t  =  .100  msec 
indicates  that  a  rather  strong  depletion  occurs  in  the  regions  closest  to 
the  walls  where  the  concentration  closely  follows  the  source's  spatial 
distribution.  This  simple  analysis  justifies  the  use  of  an  "enhanced" 

diffusion  rate  in  the  kinetics  model. 


DIFFUSION  RATES  FOR  EACH  REGION  (1mm  SPACING) 

l/r0(12  msec)  1/Tq(100  msec) 
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